Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 25 November 2009 (MN KTeX style file v2.2) 



Properties of warm absorbers in active galaxies: a systematic 
stability curve analysis 

Susmita Chakravorty^ *, Ajit K. Kembhavi^*, Martin Elvis^*, Gary Ferland^* 

^lUCAA, Post Bag 4, Ganeshkhind, Pune4U 007, India; 
^Harvard-Smithsonian Center for Astrophysics, Cambridge, MA 02138; 
^Department of Physics and Astronomy, University of Kentucky, Lexington, KY 40506 



25 November 2009 



ABSTRACT 

Signatures of warm absorbers are seen in soft X-ray spectra of about half of all Seyfertl 
galaxies observed and in some quasars and blazars. We use the thermal equilibrium curve to 
study the influence of the shape of the ionizing continuum, density and the chemical compo- 
sition of the absorbing gas on the existence and nature of the warm absorbers. We describe 
circumstances in which a stable warm absorber can exist as a multiphase medium or one with 
continuous variation in pressure. In particular we find the following results: i) the warm ab- 
sorber exists only if the spectral index of the X-ray power-law ionizing continuum a > 0.2 
and has a multiphase nature if a ~ 0.8, which interestingly is the spectral index for most of the 
observed Seyfert 1 galaxies; ii) thermal and ionization states of highly dense warm absorbers 
are sensitive to their density if the ionizing continuum is sufficiently soft, i.e. dominated by 
the ultraviolet iii) absorbing gas with super-Solar metallicity is more likely to have a multi- 
phase nature; iv) the nature of the warm absorber is significantly influenced by the absence of 
iron and associated elements which are produced in the later stages of star formation history 
in supernovae of type la. 

Key words: quasars: absorption lines - galaxies : active - Seyfert - ISM: abundances X-rays: 
ISM 



1 INTRODUCTION 

Absorption due to highly ionized Oxygen and elements of similar 
atomic number (e.g. O VII, VIII, Ne X) is commonly found in the 
soft X-ray spectra of active galactic nuclei (AGNs) (Halpern 1984; 
Nandra & Pounds, 1994; Reynolds, 1997; George etal. 1998). 
High resolution X-ray spectra from Chandra and XMM-Newton 
show that these absorption lines are always blue-shifted compared 
with the emission lines and so must be in outflow from the cen- 
tral continuum source (CoUinge etal. 2001, Kaastra etal. 2002, 
Kinkhabwala etal. 2002, Blustin etal. 2003, Krongold etal. 
2003, Netzer etaZ. 2003, Turner efaZ. 2004). The absorption edges 
and lines respond to variation of the ionizing continuum emitted by 
the central engine. For example, Krongold et al. (2007) show that 
for NGC405 1 the ionization state of the WA responds to the con- 
tinuum flux as predicted for gas that stays close to photoionization 
equilibrium, so that other factors, e.g. velocity changes and shocks, 
do not appear to be important. It is therefore reasonable to assume 
that the absorbing cloud is photoionized by the continuum radia- 
tion from the active galactic nuclei (AGN). This partially ionized 
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optically thin gas along the line of sight to the center of the AGN 
has come to be known as the warm absorber (hereafter WA). 

Understanding the structure of the WA will lead to better un- 
derstanding of the structure of AGNs; for example, identifying the 
location from which the wind arises will impose constraints on 
the physics of the wind source. Some authors find that the WA 
is well fitted with a model with two ionization phases (Krongold 
etal. 2003; Krongold etal. 2005a; Ashton etal. 2006; CoUinge 
etal. 2001) or three phases (Netzer etal. 2003; Morales, Fabian 
& Reynolds, 2000; Kaastra etal. 2002), in which the components 
are in pressure equilibrium. Other authors (Ogle etal. , 2004) and 
Steenbrugge etal. , (2005) have argued in favor of a continuous dis- 
tribution of ionization state for the WA, in objects for which others 
fit just a few phases. One way to investigate which solution is more 
likely is to investigate the physics determining the thermal structure 
of WA. 

The effect of various factors on the WA can be studied con- 
veniently using the thermal equilibrium curve (hereafter referred to 
as the stability curve) of temperature (T) against the ratio of ion- 
ization parameter (^, see below for definition) to T. Gas lying off 
the stability curve will cool or heat until reaching the curve. If the 
curve has kinks that produce multiple stable values at fixed £,/T, 
then the WA can have multiple phases in mutual equilibrium. This 
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effect may be responsible for the multiphase nature of the inter- 
stellar medium of galaxies (Field, 1965). Ionizing continuum from 
AGN can induce such kinks (Krolik, McKee & Tarter, 1981), to a 
larger extent than the interstellar medium. The pronounced kinks 
representing WA in the stability curve lie in the intermediate tem- 
perature range, 5 < log T < 7 where photoionization heating and 
collisional cooling of various metals are expected to determine the 
thermal balance of the WA. The details of the thermal stability is 
strongly affected by the shape of the ionizing continuum emitted by 
the central source (i.e. their spectral energy distributions or SEDs), 
and by the chemical composition of the absorbing gas. 

Our main motivation in this paper is to improve on previ- 
ous work (Hess, Kahn & Paerels 1997; Komossa & Mathur 2001; 
Krolik, McKee & Tarter 1981, Krolik &; Kriss 2001; Reynolds & 
Fabian 1995; Rozanska, Kowalska & Goncalves 2008) by carrying 
out a systematic investigation of the stability curves and finally to 
give a quantitative handle on the relative stability nature of the WA 
as a function of the various physical parameters involved. In Sec- 
tion 2 we discuss the physical properties of the WA. In Section 3 
we give all the constraints and assumptions used in the paper. The 
simulations of the stability curves are sensitively dependent on the 
underlying atomic physics used by the photoionization codes. We 
have shown (Chakravorty etal. , 2008, hereafter Paper 1), that di- 
electronic recombination rate coefficients have a large effect on the 
stability curve. Using recent atomic physics databases that include 
many more of these rates, with more accurate values, is important 
because it significantly changes the predictions for the physical na- 
ture of the WA. The serious impact of this effect is addressed in 
further details in Section 4 while we discuss the stability curve and 
comparison with earlier works. We choose a specific scheme to rep- 
resent the ionizing continuum (discussed in Section 3) and vary its 
parameters to create a reasonable ensemble of spectral shapes as 
demonstrated by a wide variety of AGN which harbour WA. The ef- 
fect of these different SEDs on WA nature is discussed in Section 5. 
Rozanska et al. (2008) have demonstrated that for sufficiently soft 
AGN SED, the thermal balance becomes dependent on the density 
values if it is above a certain threshold. This can only happen if 
the gas is dense enough and the ionizing continuum is appropriate 
to induce heating processes other than photoionization to dominate 
the thermal balance. Since we sample a rather large variety of ion- 
izing continua for our investigations, the present work gives us an 
appropriate opportunity to study the effect of density on stability 
curves in Section 6. Optical and ultraviolet emission and absorp- 
tion line studies suggest that central regions of AGN have solar or 
higher metallicities (Hamann & Ferland, 1999) and the WA tem- 
perature range of the stability curve is sensitive to the atomic in- 
teractions due to various metals. This abundance dependence pro- 
vides the motivation for our study of the influence of individual ele- 
ments, and groups of elements in Section 7. The improved updated 
atomic database used for the simulations provides the opportunity 
to go beyond the qualitative discussions and, for the first time, make 
quantitative comparisons of the multiphase nature and the robust- 
ness of the stability regions as a function of the ionizing SED and 
chemical composition of the absorber in Section 8. Our results are 
summarised in Section 9. 



2 WARM ABSORBERS 

Studies of WA variability in response to continuum changes show 
that it is reasonable to assume the WA to be in ionization and ther- 
mal equilibrium (Nicastro etal. , 1999). The timescales to reach 



equilibrium seem to be less than a day (Nicastro etal. 1999; Kro- 
ngold etal. 2003; Krongold etal. 2005a; Krongold etal. 2005b; 
Krongold etal. 2007). Equihbrium is reached when photoioniza- 
tion is balanced by recombination, and ionization and Compton 
heating are balanced by recombination and Compton cooling, cool- 
ing. The Compton heating and cooling become important at high 
temperature which is solely determined by the shape of the illiuni- 
nating continuum, and t5'pically is T > lO'^ K for standard Seyfert 
1 spectral energy distributions, which the clouds can attain only 
when the ionizing continuum has significant number of hard X-ray 
photons. The typical column density in WAs, estimated mainly us- 
ing absorption edges, is A^h ~ 10^^^^ cm~^ (see the references in 
Section 1). In some cases lines are easily detected whereas edges 
are not seen (Kaastra eta/. 2002; Rozanska etaL 2004). However, 
the A^H values derived from edges alone are systematically high 
by large factors (~10) - see Krongold etal. 2003, Kaspi etal. 
2002. Modeling of the absorber allows the temperature and ion- 
ization parameter to be calculated. Typically, the temperatures are 
of the order of 10'' K and the ionization parameter ^ ~ 100 (see 
definition of ^ below). However, in some objects WAs are found 
to have additional components which exhibit higher temperatures, 
~ 10® K (see the references in Section 1). 

The estimated distance of the WA from the central source cov- 
ers a wide range: from the kiloparsec scale of the extended narrow 
emission line region (Kinkhabwala etal. 2002), through distances 
of ~ 0.01 — 1 pc, i. e. , at the scale of the obscuring molecular 
torus (Krolik & Kriss, 2001; Crenshaw etal. 2003; Blustin etal. 
2005), to within the H (3 broad emission line region, on the scale 
of the accretion disk, which is a few thousand Schwartzchild radii 
(Krongold etal. , 2007), where the WA would form the base of an 
accretion disk wind (Murray & Chiang 1995, Elvis 2000). The pho- 
toionization state of the WA can be parametrised by the ionization 
parameter which is the ratio of the ionizing photon flux to the gas 
density (Tarter, Tucker & Salpeter, 1969): 

^ = Lion/nnR^ [ergcms"^], (1) 

where Lion is the luminosity between 1 — 10'^ Rydberg. Hence 
^/T ~ L/pR^, p being the gas pressure in the WA. Although, in 
this paper we have used ^ for all our discussions, in the context of 
warm absorbers, there is another popular definition of the ioniza- 
tion parameter given by 

C7 = Qion/47rcnHJi^ (2) 

Qion (s"'^) is the rate of incident photons above 1 Rydberg and c is 
the velocity of light. For a given shape of the SED of the ionizing 
continuum, the conversion from ^ to (7 is unique and can be done 
conveniently. Although luminosity can be constrained from sepa- 
rate independent observations, it is difficult to remove the degener- 
acy between density and distance in the product nn featuring in 
the denominator of the ionization parameter. Thus, the estimation 
of distances for WA becomes challenging; for example, Netzer et 
al. 2003 and Krongold et al. 2006 put upper limits differing by up 
to three orders of magnitude on the WA distance in the same object 
NGC3783 (but using different observations). As the mass loss rate 
is proportional to the radius such large factors of uncertainty have 
a large effect on AGN energetics and on feedback from the AGN to 
the surrounding media. 

Krohk, McKee and Tarter (1981) showed that a stable state 
exists at T 10* K for a photoionized cloud which undergoes 
collisional cooling. There can be an additional stable state for the 
absorbing cloud at T ~ 10® K, if the incident photoionizing con- 
tinuum has a hard X-ray component as an extra heating agent. Thus 
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they suggest that presence of a hard X-ray continuum would result 
in clumping, with cold clouds remaining in equilibrium with a hot 
medium. However, for optically thin, ionized gas, the temperature 
range 10^ — 10^ K is well known to be an imstable zone, bringing 
into question the model of WA as gas in thermal equilibrium. Ra- 
diation in spectral lines is the primary cooling mechanism in this 
temperature range. The contribution from radiative recombination 
and two-photon continuum emission (Gehrels & Williams 1993) 
form some additional small regions of stability. 



3 CONSTRAINTS AND ASSUMPTIONS 
3.1 The code and its implementation 

We assume a simple model of a geometrically thick, optically thin, 
plane parallel slab of gas illuminated on one face for the absorb- 
ing medium which is located at a certain distance R from a cen- 
tral continuum source of bolometric luminosity L. The geometrical 
thickness of the gas is taken to be much less than R. The physi- 
cal condition of the absorbing medium is determined by the flux 
received at the illuminated face of the gas {L/A-kR?), the shape 
of the ionizing continuum (SED), the number density of the gas 
riH, the column density of the gas A^h, and the chemical compo- 
sition of the absorbing gas. To obtain the equilibrium conditions 
for a gas with an assumed set of parameters which specify these 
inputs, we have used the publicly available photoionization code 
cLouDY^version C07.02 (hereafter C07), see Ferland (1998). 

CLOUDY calculates the ionization equilibrium conditions by solv- 
ing the energy and charge conservation equations under the as- 
sumption that all the atomic processes have had time to reach a 
steady state. The density rii of i*** species is then given by the bal- 
ance equation 

dri; 
'di 



= n-jRji + Source 



(3) 




+ Sink 



0, 



where Rji is the rate of species j going to i. Source is the rate per 
unit volume of appearance of new atoms in i*'' species and Sink is 
the rate that they are lost (Osterbrock & Ferland, 2006). Thermal 
equilibrium conditions are calculated by evaluating the heating and 
cooling rates at different temperatures for the given gas and finding 
the temperatures where the two rates balance each other. 

For all the calculations in this paper we have used cloudy in the 
constant density mode which is chosen for the sake of numerical 
simplicity, as we probe a rather large range of parameter space to 
determine the properties of the WA. We would like to note here that 
using the TITAN code (Dumont, Abrassart & CoUin, 2000), Rozan- 
ska etal. (2006) and Goncalves etal. (2006) have pioneered the 
constant total pressure mode of calculations for photoionization 
equilibrium where they relax the approximation that the density is 
constant. Their only assumption is that the line of sight gas is in 
total pressure equilibrium and they introduce additional iterations 
between temperature and density within the code in a self consis- 
tent manner. This results in natural stratification of the cloud with 
different temperatures and densities where each layer is illuminated 



by the radiation already absorbed by the previous layer. These two 
methods agree well with each other for optically thin clouds, but 
show significant differences when the column density of the ab- 
sorbing cloud is high (iVn > 10^^). In our future publications 
where we will restrict ourselves to study equilibrium properties of 
WA as functions of a smaller range of physical parameters, we will 
use similar methods. 



3.2 Density dependence 

Estimation of temperature is independent of nn, for a gas in ther- 
mal and ionization equilibrium when the dominant processes of 
heating and cooling are photoionization and recombination respec- 
tively. This is because in this regime of thermal equilibrium, den- 
sity enters the temperature estimation only through the ionization 
parameter. Hence the stability curve would be insensitive to the 
value of the density used. This gives us the liberty to choose an 
intermediate value for WA density, nu ^ l(f cin^'\ loading to 
L/R^ ~ 10^^ from Equation 1, for typical values ^ ^ 100 corre- 
sponding to WA. Except for Section 6, where we check the robust- 
ness of this assimiption, we have assumed the value lO'^'^ for L/R^ 
throughout the paper. 



3.3 Ionizing continuum 

The 1 — 10 keV part of the observed spectra of AGN is usually 
modeled as a simple power-law with exponential lower and upper 
cut-offs specified by fmin and Umax respectively. 



(4) 



where /{v) is the flux per unit frequency interval and a{> 0) is the 
energy spectral index. For most quasars, the X-ray spectral index 
lies in the range 0.7 < a < 0.9 as shown by Wilkes & Elvis in 
1987 and also recently by Grupe etaZ. 2006 and Lopez etaZ. 2006. 
However some objects show extreme values like a < or a > 2. 
There is a systematic difference between radio-quiet and radio-loud 
AGNs in that radio-quite AGNs have steeper X-ray spectra (a ~ 
1.0) than radio-loud ones, (a ~ 0.5, Wilkes & Elvis 1987; Grupe 
etal. 2006; Lopez etal. 2006). 

Ionizing continua given by Equation 4 have been used earlier 
by Reynolds & Fabian (1995, hereafter RF95) for doing detailed 
analysis of nature of WA as a function of various relevant physi- 
cal parameters. Their work is further discussed in Section 4. The 
dashed and dotted curve labeled as Early set in the lower panel of 
Figure 1 is drawn using Equation 4 with oc = 0.8, Emin = 13.6 eV 
and Emax = 40 keV. 

The extension of the soft X-ray SED into the extreme ultravio- 
let (EUV) cannot be directly observed in most cases because energy 
in EUV is easily absorbed by our own Galaxy. The join between the 
EUV and X-ray flux is parametrized by the slope aox defined as 



-0.384 log 



F^{2keV) 



F^{2500A) 



(5) 



(Tananbaum etal. , 1979). To obtain aox in the observed range 
1 < aox < 2 we need a steeper soft component in addition to the 
power-law given by Equation 4. We assume a two component SED 
to simultaneously include the information on X-ray and EUV and 
which can be expressed as 
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for z/ > 2500A, 



(6) 



© 0000 RAS, MNRAS 000, 000-000 



4 Chakravorty et al. 



where as(> 0) is the spectral index of the steep soft component 
and r) is the relative normalization factor. For energies lower than 
2500A we have used a spectral shape given by 

f{u) ~ i/""-^ for 0.206 Ryd <v < 2500A 

~ i/"^" for 9.f2 X 10"^ <v < 0.206 Ryd 
~ v^-'" for u < 9.f2 X fO~^ 



(7) 



which is similar to the shape used by Mathew & Ferland (1987) 
for the same energies. Equation 6 and Equation 7 taken together 
describe a more realistic ionizing continuum which irradiates the 
WA. Such a continuum is shown as the solid curve labeled Standard 
set in Figure 1 with a = 0.8, Qs = 2.0, aox ~ 1-2, Emax = 
200 keV and lower energy cut-off following Equation 7. 

In the lower panel of Figure 1, the Standard set continuum is 
compared with various other SEDs used or discussed in this paper 
and we will refer to these comparisons in more details in Section 4. 
In this paper we have used the simple definition for the ionizing 
continuum, as in Equation 4, only in Section 4 where we have 
given an introductory description of the stability curve and com- 
pared results with similar works in literature (RF95). In all other 
sections we have considered the more realistic Standard set spec- 
tra. All through the paper we have assumed Os = 2, and for a given 
value of Q, obtained the desired value of aox by adjusting rj. In Sec- 
tion 5, when we investigate the changes in the nature of the stability 
curve with variation in the shape of the continuum, we consider a 
wide range of slopes, 0.2 < a < 1.1 and 1 < aox < 1-8 and a 
wide range of higher energy cutoff, 50 keV < Emax < 400 keV 
to incorporate various types of AGN. 



4 STABILITY CURVE AND EARLIER STUDIES ON IT 
4.1 Definition 

The stability curve, may be considered to be a phase diagram in the 
log T — log(5/r) plane, with each point on the curve representing 
an equilibrium state of the system at a temperature T, for ioniza- 
tion parameter ^. The other parameters describing the system have 
specific values as mentioned below. The solid curve in Figure 2 
is generated using version C07 of cloudy, for an absorber having 
gas density nu = 10^^/^ and solar metallicity, being ionized by 
a power-law continuum having energy spectral index a = 0.8 and 
extending from 13.6 eV to 40 keV; hereafter we will refer to this 
as the Early model and the corresponding parameter set as the Early 
set. Computations were done using a grid of models with a range of 
values of the ionization parameter 5, and the code is constrained to 
perform single zone calculations to ensure that the gas is optically 
thin. The ionizing continuum used in the Early model is shown as 
the dotted-and-dashed curve in the lower panel of Figure 1 . We use 
the Early set in this section to compare results with earlier works 
like RF95 who have done extensive study of the stability properties 
of the WA as a function of various shapes of the ionizing contin- 
uum. 

An isobaric perturbation of a system in equilibrium is repre- 
sented by a small vertical displacement from the stability curve; 
such perturbations leave C/^ ~ L/pB? constant, which for con- 
stant L/Ft^ leave the pressure unchanged. If the system is located 
on a part of the curve with positive slope, then a perturbation cor- 
responding to an increase in temperature leads to cooling, while a 
decrease in temperature leads to heating of the gas. A gas in such a 
state is therefore thermally stable. On the other hand, a gas located 
in a region of the curve with negative slope is thermally unstable, 
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Figure 1. Top panel : Ionizing continuum obtained using Equation 6 and 
Equation 7 taken together, where a = 0.8, «a = 2.0, aox = 1-2 and the 
continuum is exponentially cut-off at 200 keV. At lower energies, this con- 
tinuum is cut-off following Equation 7 with the successive spectral indices 
of —0.5, —1.0 and 2.5. This continuum is labeled as the Standard set. We 
have also appropriately labeled the dominant component of the SED at vai- 
ious energies and used the dashed line to show the resultant aox = 1-2. 
Bottom panel : Comparison of the Standard .set ionizing continuum with 
vaiious other SEDs used or discussed in this paper The dotted-and-dashed 
line labeled as the Early set is the ionizing continuum given by Equation 4 
with a X-ray slope a = 0.8 and exponential cut-offs at Emin = 13.6 eV 
and Emax = 40 keV. The dotted line represents the ionizing continuum 
used by Krolik etal. 1981 and the dashed Hne that used by Mathew & 
Ferland (1987). The important energy values including the upper and lower 
energy cut-offs (13.6 cV, 40 keV and 200 koV) and the range of definition 
for aox (2500A and 2 keV) have been marked and labeled. 



and may achieve a multiphase equilibrium under certain physical 
conditions which we will discuss in later sections. 



4.2 Earlier studies 

Before we go on to discuss our own results with the stability curve, 
here, we would like to mention some of the earlier works done by 
various authors on WA and/or stability curves. 

Krolik etal. (1981) obtained the stability curve for AGN 
using parameters similar to our Early set except for the ionizing 



© 0000 RAS, MNRAS 000, 000^00 



Stability curves for Warm Absorbers in AGN 5 



continuum. The continuum used by them is plotted as the dotted 
curve in Figure 1. They have argued in favour of ~ 10^ K stable 
states remaining in pressure equiUbrium with the ~ 10* K phases. 
Although there are some narrow stable phases at ~ 10^ K, they 
have ignored the significance of such states. Moreover their stabil- 
ity curves do not show any stable states at ~ 10® K. 

In a similar study on the absorbing medium, Mathew & Fer- 
land (1987) do not find the ~ 10"* K and ~ 10* K phases in pres- 
sure equilibrium with each other. The continuum used by Mathew 
&. Ferland (1987) is widely used as a broad band continuum for 
AGN modeling studies and has been represented with the dashed 
line in Figure 1 . The energy spectral index of the ionizing contin- 
uum in the Standard set is very similar to that of the Mathew & Fer- 
land SED in X-rays (0.3 - 100 keV). However, it is not objective 
to make a comparison at this stage because, although the ionizing 
continuum in the Standard set given by Equation 6 represents the 
overall shape of the SED between 3.6 eV to 2 keV using the def- 
inition of Qoa;, it is unable to give information on any underlying 
specific features in this energy range. In this paper we will constrain 
ourselves to use ionizing continuum given by Equation 6 and Equa- 
tion 7, but will investigate the details of the spectral shape in the en- 
ergy range (I - 1000 eV) in subsequent publications (Chakravorty 
etal. , in preparation) by incorporating the effect of soft excess at 
~ 150 eV and a disk blackbody component at ~ 20 eV. 

Reynolds and Fabian (1995, RF95) have studied the warm ab- 
sorber stability conditions using physical conditions as described 
by our Early set of parameters. They used a simple power-law for 
the ionizing continuum and varied the power-law index ot to check 
for its effects on the WA and also examined some basic effects of 
adding a soft excess component to the power-law ionizing contin- 
uum. Using power-law continuum with a — 0.8, they find that 
the stability curve is independent of density variation in the range 
10 < < 10" cm"^ 

Hess et al. (1997) address the discrepancy that the physical 
properties inferred from observations of some low-mass X-ray bi- 
naries and Seyfert galaxies correspond to thermally unstable re- 
gions of the stability curve predicted by photoionization codes. 
They have carried out an extensive study of the stability curve as a 
function of various ionizing continuum and different abundance of 
the absorbing gas. They have studied the cause of the instabilities 
in great detail and have cautioned against the outdated fits and pa- 
rameters used by the simulation code and also against the accuracy 
of the atomic data then used for determining the thermal balance. 

Komossa & Meerschweinchen (2000) and Komossa & Mathur 
(2001) have extensively studied the stability curve respectively as 
a function of ionizing continuum and chemical composition of the 
absorber with the aim of explaining the WA in various narrow line 
Seyfert I galaxies. They find that steeper ionizing continuum sta- 
bilizes WA and super-Solar metallicity of the gas enhances the 
possibility of multiphase WA. By modeling the observations of 
NGC4051, indications are found that the WA in this object has 
super-Solar abundance. 

4.3 Comparison with earlier work: importance of updated 
atomic data 

In the WA temperature range 10^ < log T < 10^ K cooling pro- 
cesses are dominated by coUisional recombination which can be 
radiative or dielectronic. These processes lead to local regions of 
thermal stability in this otherwise unstable temperature range as 
discussed by Gehrels & Williams (1993), Hess, Kahn & Paerels 
(1997) and also depicted by RF95. We can see these effects in our 
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Figure 2. Stability curves showing the distribution of equilibrium temper- 
ature T as a function of ^/T. It is assumed that an optically thin shell of 
gas of solar metallicity is ionized by a continuum with energy a = 0.8, ex- 
tending from 13.6 eV to 40 koV. The regions of the plane where heating or 
cooling dominates are indicated. The solid curve is obtained using version 
C07 of CLOUDY while the dotted curve is for version c8412a. The Early set 
of parameters, as described in the text, is used in generating the two curves, 
which are seen to be significantly different for lO*'® < logT < 10^ K. 



stability curve generated using the Early set of parameter values 
and shown as the solid curve in Figure 2. However the atomic 
physics data base has been updated considerably over the last 
decade and our results show significant quantitative variations from 
earlier work which we emphasise by making explicit comparisons 
with RF95. 

The dotted line in Figure 2 is a reproduction of the stability 
curve (using version C8412a of cloudy) in Figure 3 of RF95 where 
they had used the Early set of physical parameters. The variations 
in the C07 stability curve from the C84 curve lead to different phys- 
ical predictions for the WA; for example, the 10'^ K state becomes 
much more stable in our work resulting in more pronounced pos- 
sibility of multiphase nature of the WA. The 5 < logT < 6.5 
region of the stability curve, relevant for WA, sensitively depends 
on the detailed atomic physics of the various elements which con- 
tribute to photoelectric heating and cooling due to recombination 
and collisionally excited lines. In Paper I we have shown that the 
total recombination rates (dielectronic -I- radiative), for the signif- 
icant cooling agents which bring about the maximum variation in 
the stability curves in Figure 2, are predicted to be larger in C07 
than in C84. 

In the WA temperature range lO'"' < T < lO'^ K, dielectronic 
recombination dominates over radiative recombination for many 
ions (Osterbrock & Ferland, 2006). Unlike the radiative recombi- 
nation rate coefficients, dielectronic recombination rate coefficients 
(DRRC) have undergone significant changes over the last decade. 
The importance of revisiting the values of DRRC were first pointed 
out by Savin et al. (1997) and further emphasised by a series of pa- 
pers by Badnell and coworkers (Colgan, Pindzola & Badnell 2004, 
Colgan et al. 2003, Altun et al. 2004, Zatsarinny et al. 2004a, Mit- 
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Model 



Definition 



Ionizing continuum 



Low energy 
Cutoff 



riH 



MetaUicity 



Early 



Equation 4 



0.8 



40keV 



(Equation 4) 



10^ cm-3 



Zone 1 
(See text) 



Standard 



Equation 6 
Equation 7 



0.8 



1.2 



2.0 



200 keV 



Equation 7 



io"/€ 

(See text) 



10^ 



Table 1. Comparison of the various parameters used to constitute the Early and the Standard model for stabUity ciu^es. 



nik & Badnell 2004, Zatsarinny et al. 2003, Gu 2003, Zatsarinny et 
al. 2004b and Gu 2004). The DRRC data base is more extensive 
now and shows that the rate coefficients for different ionization 
states of various elements are substantially larger than the values 
used previously, but the database is still not complete, especially 
for the lower ionization states. For ions which do not have com- 
puted DRRC values, C07 uses a solution, as suggested by Ali et al. 
(1991); for any given Icinetic temperature, ions that lacis data are 
given DRRC values that are the averages of all ions with the same 
charge. The advantage of this method is that the assumed rates are 
within the range of existing published rates at the given kinetic tem- 
perature and hence cannot be drastically wrong. However, in Paper 
I, we have checked whether the points in the C07 stability curves 
have computed DRRC values or guessed average values, concen- 
trating on the parts of the curve which have multiphase solutions 
for the warm absorber (10°^ < T < lO" K) and are different from 
the C84 stability curve. We find that all the ions which act as major 
cooling agents for each of these points in the stability curve have 
reliable computed DRRC values which are not likely to change in 
near future. Thus, the new data base provides a more robust mea- 
surement of the various physical parameters involved in studying 
the thermal and ionization equilibrium of photoionised gas. This 
not only gives us the motivation to improve upon previous qual- 
itative stability curve analysis by various authors (Hess, Kahn & 
Paerels 1997; Komossa & Mathur 2001; Krolik, McKee & Tarter 
1981; Krolik &; Kriss 2001; Reynolds & Fabian 1995; Rozanska 
et al. 2008) but also provides the chance to extend their work with 
quantitative estimates on the nature of the multiphase WA in Sec- 
tion 8. 



4.4 Standard set 

The ionizing continuum used in the Early set has now been shown 
by improved spectral UV and X-ray measurements not to corre- 
spond to observed continua (Zheng etai. 1997, Shang etaZ. 2005). 
As discussed in Section 3, we will use a more realistic SED with 
two power-law components given by Equation 6 and lower en- 
ergy cut-off given by Equation 7. The fiducial value for aox for 
Seyfertl galaxies is 1.2 (Netzer 1993). We define a Standard set of 
parameter values to be used in the rest of the paper where the ion- 
izing continuum shown as the solid line in Figure 1 is parametrized 
by a = 0.8, as = 2 and r? is appropriately adjusted to attain 



aox = 1.2. Emax for both components is 200 keV. The high en- 
ergy cut-off is an average of the observational values reported from 
Beppo-Sax which found hard X-ray spectral cut-off to span the 
range 50 to 450 keV (Matt 2004; Tueller et al. 2008). The absorb- 
ing gas in this Standard model is assumed to have solar abundance 
{Z — Zq) as given by Allende Prieto, Lambert & Asplund (2002, 
2001), Holweger (2001) and Grevesse & Sauval (1998). The den- 
sity is given by nn = 10^^ as in the Early set. In addition we 
assume A'h = 10^^ cm~^, for the slab of the gas in the Standard 
set. Table 1 shows a comparison of the various parameters used in 
the Early set and the Standard set. 



5 INFLUENCE OF EUV/ SOFT X-RAY SLOPE 

In this section we consider the effect on the stability curve of chang- 
ing the shape of the continuum in Equation 6 by varying a and aox 
over a range of values. All through the section we assume that the 
gas has chemical composition of solar abtmdance. 



5.1 Influence of Soft X-ray slope 

The left panels of Figure 3 shows the shapes of the ionizing con- 
tinua and the corresponding stability curves for different values of 
a, aox is held at a constant value of 1.2 for these curves. We have 
highlighted the stable 10'''' K < T < 10^ K warm absorber states. 

For very flat continua, with energy index a ^ 0.2 or less, al- 
most no stable state exists in the temperature range of 10* — 10^ K, 
since there are no regions of the curve here with positive slope ex- 
cept for a very narrow one at log(^/T) ~ —2.6. The large num- 
ber of high energy photons pushes the Comptonization tempera- 
ture to T > 10* K. As noted by Krolik etal. (1981), there is a 
cooler phase at log T < 4.5 and a hot phase which is transparent 
to atomic features. Warm absorbers are thus unlikely to be found 
in AGN which have such flat X-ray spectra. Stability curves for 
steeper incident spectra show additional stable states in the inter- 
mediate temperature regime. Curves corresponding to a = 0.5 and 
a = 0.8 are both stable up to ~ 2.5 x 10® K. The standard curve 
shows one extensive stable state at ~ 10® K. We will discuss the 
effect of multiple phases further in Section 8. For steeper spectra 
with a = 1.1, the absorbing gas is predicted to be stable for all 
^ with the Compton temperature now lying at log T = 6.9. The 
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Figure 3. The shapes of the ionizing continua and the corresponding stability curves for varying X-ray slope a (left) and valuing EUV to X-ray slope, oiox 
(right) for an optically thin shell of gas of solar metallicity. In the top panels showing the SEDs, the important energy values including the upper energy cut-off 
(200 keV) and the range of definition for (2500A and 2 keV) have been marked and labeled. Left panels: X-ray flux index a is varied from one curve to 
the other in the range 0.2 < o < 1.1 with aox held constant at the value 1.2. Right panels: The EUV to soft X-ray slope Oox is varied over a range of values 
with the X-ray slope a held constant at the standard value 0.8. On the stability curves we have used black on green, blue on red, gray on black, red on blue 
and cyan on magenta to highlight the stable WA phases. Note that in the left panel, the log T axis extends up to 9. The flat spectrum (a = 0.2) curves have 
higher temperature Compton stabihzation at T > 10* K. 



distinct multiple phases at different temperatures, but in near pres- 
sure equilibrium are lost. For such steep continua, stable states exist 
over a continuous distribution of ^/T and T in the phase space. 

In summary, AGN having flat X-ray spectra with a < 0.4 are 
not likely to have WAs. Although a steeper continuum can lead to 
physical conditions that can sustain warm absorbing gas in thermal 
equilibrium, multiphase nature of WA can only be seen for AGN 
spectra with intermediate slopes of a ~ 0.8. It is interesting to note 
that multiphase WAs are predicted for continua having slopes in the 



range 0.7 < ct < 0.9 which is the observed range for the bulk of 
the observed AGN population (Wilkes & Elvis 1987; Grupe etal. 
2006; Lopez efaL 2006). 

We would like to note an interesting difference between re- 
sults obtained by us and that by Ki^olik & Kriss (2001, hereafter 
KKOl). Using a model ionizing continuum given by Kaspi et al. 
(2001), KKOl have obtained stability curves computed by version 
2.1 of XSTAR (Kallman 2000^^). The ionizing continuum by Kaspi 
et al. (2001) has values of oiox ~ 1-3 and a = 0.77, very sim- 
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ilar to those used in the Standard set. However, KKOl find that 
the corresponding stabihty curve shows a steep near vertical rise in 
temperature from ~ 3 x lO" K to ~ lO'' K for 40 < C < 1500 
which gives a continuous distribution of stable states over a wide 
range of ^, while pressure (^/T) remains almost constant. For the 
same range of (50 — 1200), our Standard curve, shows two dis- 
tinct WA phases in pressure equilibrium with each other and the 
change in temperature is gradual and spans a much wider range of 
C/T values. 

5.2 Influence of aox 

In the right panels of Figure 3 we have shown the results of vary- 
ing the EUV to soft X-ray slope aox (see Equation 5). The X-ray 
slope a is held at a constant value of 0.8. Observations find aox be- 
tween ^ 1 and 2 as reported by Strateva etal. 2005 and Steffen 
etal. 2006. They find that there is an anti correlation between aox 
and the ultraviolet luminosity Lj/(2500A) as was first reported by 
Tananbaum et al. in 1979. We have used the representative values 
of 1.0, 1.4, 1.6 and 1.8 in addition to the Standard value 1.2. The 
WA phases are highlighted as before. 

From Figure 3 we see that the nature of stability in the ctox ~ 
1 curve is very similar to that in the aox = 1-2 curve. However, 
beyond 1.2, with a sharper drop in flux from the EUV to the soft 
X-ray region, i. e. , with steeper value of Oox, there is a gradual in- 
crease in the stability of the absorbing gas at 10^ K, and a gradual 
shift to higher pressure (^/T). The aox = 1.4 curve shows three 
distinct stable states at 10^ , 10® and 10^ K in pressure equilibrium 
with each other. For steeper Oox = 1-6, 1.8, stability at 10^ K is in- 
creased. The aox = 1.8 curve shows no unstable states. For such 
spectra, the WA can have a continuous distribution of pressure and 
density. The Compton temperature for the stability curves under- 
standably gradually drops from logT ~ 7.7 for aox ~ 1.0 to 
log T ~ 6.0 for aox ~ 1.8 as the ratio of EUV to X-ray photons 
increases and Comptonization is progressively dominated by lower 
energy photons. 

To summarise, the existence of the WA is independent of the 
EUV to X-ray slope of the ionizing continuum. However, its multi- 
phase nature is influenced by this slope; the steeper the fall of flux 
from EUV to X-ray (i. e. , larger aox, the lesser is the possibility of 
finding a multiphase WA. 

5.3 Influence of Emax 

The Standard set ionizing continuum is exponentially cut off at 
Emax = 200 keV (see Equation 6). Results from the Swift/BAT 
(Burst Alert Telescope) hard X-ray sky survey (Tueller et al., 2008) 
show that, of the brightest few dozen AGNs, for which Emax can 
be determined, the lowest value is 40 keV, while some 2/3 of the 
objects have Emax > 120 keV, the highest value that BAT can 
measure. As a result an appropriate value of Emax for the Stan- 
dard Set needs to be set well above the Early Set value. We adopt 
Emax ~ 200 keV. The high energy cut-off is important in Comp- 
ton processes and hence in determining the Compton temperature 
of a stability curve, but is not likely to affect its WA region. How- 
ever, for the sake of completeness we examine the effect on stability 
curves of varying the value of Emax- 

Figure 4 shows the various stability curves drawn using the 

^ URL : http://heasarc.gsfc.nasa.gov/docs/software/xstar/docs/html/ 
xstarmanual.html. 




-4 -3.5 -3 -2.5 -2 

log(^/T) 

Figure 4. Stability curves drawn using different values for the exponential 
higher energy cut-off Emax ■ The curves ai'e appropriately labeled. All other 
physical parameters are same as in the Standard set 



different values for Emax -The curves are appropriately labeled. All 
other physical parameters are as in the Standard set. For the lowest 
considered value of Emax = 50 keV, the Compton processes start 
dominating the thermal balance at lower temperatures resulting in 
enhanced stability at T ~ 10® K. However, the influence is still 
not significant enough to change the overall qualitative nature of 
the WA or its multiphase properties. The majority of the Swift/BAT 
AGN have Emax > 100 keV. For such values of Emax, the WA 
remains effectively insensitive to the variation in the high energy 
cut-off; the only influence is in increasing the Compton temperature 
with the increase in Emax. 



6 INFLUENCE OF DENSITY OF THE ABSORBER 

As stated in Section 3, we have assumed that the stability curve is 
insensitive to the variation in the value of hydrogen density nu of 
the WA gas. Based on this assumption, we have taken the liberty to 
choose an arbitrary fiducial value of nu ~ lO'"* cm^'^ for ^ = 100, 
the typical value of ionization parameter reported for WAs, thus 
yielding L/R^ — 10^^ from Equation 1, which we hold constant 
throughout this paper except in this section. Many authors carrying 
out photoionization calculations make similar assumptions (Hess, 
Kahn & Paerels, 1997; Krolik, McKee & Tarter, 1981; Krolik & 
Kriss, 2001, RF95) and RF95 have even verified this statement over 
a wide range of densities 10 < wh < 10^^ cm~^ for a simple 
power-law ionizing continuum with a = 0.8. However, more re- 
cently, Rozanska, Kowalska & Goncalves (2008) showed that the 
stability curve shows remarkable density dependence if the ionizing 
continuum has a complicated shape with relatively more soft X-ray, 
photons although the same exercise with a simple power-law ion- 
izing continuum with a = 1.0 yields no density dependence on 
the part of the stability curve. In this paper we have studied the sta- 
bility curve's behaviour as a function of complicated ionizing con- 
tinua with having different ratios of the soft to hard X-ray photons. 
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Figure 5. The density dependence of the stabiHty curves as a function of the ionizing continuum. The top panels correspond to the various values of the X-ray 
slope a whereas the bottom panels are for different values of EUV to X-ray slope apx and all the panels have been appropriately labeled. Density at any 
point of a stability curve is given by uh = {L / B?) / For each ionizing continuum, several stability curves are drawn by varying the ratio L/B? from 
10^ to 10^5 in multiplicative steps of 10 such that the density at ^ = 100 varies from nj] = 10^ to 10^"^ cm^'^. The dotted black lines cutting across the 
curves are the constant ^ = 100 lines. In the range nn = 10^ to lO'^'' cm""^, the density value does not affect the WA range of the stability curve at all. For 
uh > 10^" cm~'^, the density dependence of the stability curve is seen to increase with the steepness of the ionizing continuum. 



Hence it is an interesting exercise to check if the stabihty curves, 
especially in the range lO'' < T < 10^ K, drawn using the differ- 
ent ionizing continua considered in this paper show any variation if 
the density is changed. 

In this section we relax the assumption made in Section 3 
that L/R^ — 10^^ and for each ionizing continuum considered 
in Section 5, draw several stability curves by varying this ratio 
from 10^ to 10^^ such that the density at ^ = 100 varies from 
riH = 10^ to 10^^ cm^"^ in multiplicative steps of ten from one 
curve to another. The results are shown in Figure 5. The upper pan- 
els in the figure show the stability curves for continua with dif- 
ferent values of the X-ray slope a = 0.5, 0.8 and 1.1, whereas 
the bottom panels correspond to different values of the EUV to X- 
ray slope aox = 1.4, 1.6 and 1.8. The determination of thermal 
equilibrium of the WA is seen to be independent of density in the 
range 10^ < nn < lO^^cm"^, whatever be the shape of the ioniz- 
ing continuum. Stability curves with tih = 10^^ and 10^^ cm^"^ 
show insignificantly small differences from the rest of the curves 
with lower densities as we vary ol from 0.5 through 0.8 to 1.1 This 
indicates that the soft X-ray slope has little or no role to play in 



inducing density dependence to the estimations of thermal equilib- 
rium conditions. The most interesting result of this exercise is that 
as the EUV to X-ray slope oiox is increased from 1.2 through 1.4 
and 1.6 to 1.8, the WA thermal equilibrium properties become sig- 
nificantly density dependent for riH > 10^^ cm^^ and the depen- 
dence is more pronounced with the increase in olqx i.e. with the 
increase in the relative number of EUV photons, ^ 10^^ cm""^ 
can be considered to be an upper limit on WA density and would 
be consistent with clouds residing within the H /3 broad emission 
line region, on the scale of the accretion disc, a few thousand 
Schwartzchild radii (Krongold etal. , 2007). Thus higher values 
of density are unphysical for WA. However, in each of the panels 
in Figure 5, we have also plotted the nn ~ 10^^ cm"'^ curve for 
comparison and it shows gross difference from the stability curves 
with lower densities for all the continuum shapes considered. 

Since the olqx = 1-8 curves show the maximum sensitivity 
to density variation, we select them to investigate the cause of this 
effect. We choose the ^ = 100 model from each of the curves 
with densities 10^ < nn < lO'^^ cm~^ and check what are the 
dominant heating and cooling agents in determining the thermal 
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balance. The fractional contribution towards the total cooling rate 
by the dominant cooling agents undergoes only a gradual change 
as we change the density though the above mentioned range. The 
remarkable change is shown by the dominant heating agent. For 
uh < 10^" cm^'', the heating agents are all photoionized species 
and free-free heating starts to dominate for nn ~ 10^^ cm""^. The 
percentage of heating by free-free absorption increases from ~ 17 
for nn = lO" cm"^ to - 83 for uh = 10^^ cm"^. 

The ionization equilibrium conditions of the absorbing gas can 
be uniquely determined if the shape of the ionizing continuum and 
the flux received by the gas is specified. Assuming that the absorber 
is in ionization equilibrium, cloudy evaluates the heating and cool- 
ing rates at different temperatures for the given gas and determines 
the thermal equilibrium solution for the temperature where these 
two rates balance each other. In the WA temperature regime where 
photoionization is the dominant heating process and thermal equi- 
librium is achieved by recombination cooling, the temperature es- 
timates become independent of density as both the heating and the 
cooling terms have similar dependence on density which cancels 
out from both sides of the thermal balance equation. At higher tem- 
peratures, the thermal balance is attained by Compton heating and 
cooling, in which case also, the density cancels out. 

However, the interplay between various heating and cooling 
processes is also a function of the ionization continuum. For spec- 
tral shapes with steep aox there maybe sufficiently large number 
of low energy photons to make the free-free absorption by ionized 
hydrogen dominate over photoionization as the heating process. 
The heating rate for this interaction has a different dependence on 
density and will result in making the thermal equilibrium condi- 
tions sensitive to density variations for nn greater than a certain 
threshold value. 

The interesting point to note at the end of this exercise is that 
although the mode of photoionization calculations used by us and 
Rozanska et al. (2008) are different, our results are qualitatively 
consistent with each other. For sufficiently soft UV dominated ion- 
izing continuum the thermal equilibrium calculations become den- 
sity dependent if the gas is denser than a certain threshold. Since the 
gas is also assumed to be in ionization equilibrium, the ionization 
structure will also be influenced. The effects would be much more 
pronounced for constant pressure mode of calculations as shown 
by Rozanska et al. (2008). The importance of this result is that for 
dense gas illuminated by UV dominated ionizing continua, the den- 
sity of the gas can be estimated directly by modeling the observed 
spectra in UV. 



7 INFLUENCE OF THE CHEMICAL COMPOSITION OF 
THE ABSORBER 

The behavior of the models along the stability curve, from log T ~ 
4.0 to logT ~ 7.0 is sensitive to the role of atomic physics in 
determining the thermal state of the gas. The photoionization cross 
section a^(X"'"') from the ground level of X"'"' ion with the thresh- 
old is given by the equation 



-20 



n{X+')r{X+') = n{X+') 



n{X^ 



h V 



a^(X^')dv 



R4X+')du. 



(8) 



I 



-22 - 



-24 - 



0-26 - 



-28 




2.5 



log[Energy(eV)] 



cv 






: : NeX 


CVI : 


OVII 


OVIII 


■ FeXVII 



where n{X'^') is the number density of the ion, r(X+') is 
the number of excitations per particle to the (i + l)"* state and Ji, 



Figure 6. The radiation field interaction rate is plotted as a function of en- 
ergy for physical models with different ionization parameters chosen from 
different regions of the standard stability curve. See text for details. The dot- 
ted vertical lines mai'k the absorption edges of some of the ions important 
for WA, viz. CV, CVI, OVII, OVIII, FeXVII, NeX. 



is the mean intensity. Thus the function Rv{X^^) is a measure of 
the interaction between the radiation field and the X^' ion. 
In Figure 6 we have shown the dependence of 

Rv — ^ ^ Ri, (X"*"') 

i 

on ^ and hence on the ionization state of the absorbing gas. We 
have selected three different models along the Standard stability 
curve and plotted Rv/nn as a function of energy with the values 
of log^ and logT indicated for each curve; nn = 10"/5 being 
the hydrogen density. The absorption edges of various important 
ions including CV, CVI, OVII, OVIII, FeXVII and NeX which are 
the signatures of WA have been indicated with dotted lines. As the 
gas becomes more ionized by factors of ~ 10 its interaction with 
radiation decreases by orders of magnitude, especially at energies 
less than 1 keV. At any given energy the value of R^ is lower for 
higher ^. However, with highly ionized gas, the interaction shifts to 
the higher energies. For log ^ = 1, there is no sharp increase in R^ 
at energies corresponding to the WA edges; these interactions grow 
with 5 as seen in the middle and the bottom curves corresponding 
to log ^ = 2.4 and 3.4 respectively. 

We see above that at values of log ^ relevant for WAs, the na- 
ture of the interaction between gas and radiation is significantly 
influenced by the atomic physics of the heavier elements. In the 
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Figure 7. Stability curves for absorbers with different abundances. The left panel shows curves having abundance from Zq to 9Zq in steps of 2. The curves 
show that absorbers with super-Solar chemical compositions will enhance multiphase nature of the absorber. Except for the curves in the inset window in 
the left panel, the Standard set ionizing continuum with a = 0.8 has been used for the calculations. The stabihty curves in the inset window are, however, 
obtained using an ionizing continuum with flat X-ray slope of a = 0.2 to investigate if enhanced metallicity helps sustain WA for gas irradiated with such flat 
continuum. The results show no significant stable states in the temperature range 5 < log T < 7 for WAs. The right panel shows the curves having sub-Solar 
abundances from Q.OZq to Zq in steps of 0.2. Decreases in metallicity are seen to stabilize the curve in the temperature range of lVr> - 10^ K. However, 
the possibility of multiphase media decreases with the decrease in the abundance. In all the panels the stable WA are highlighted as before. The dotted black 
lines cutting across the stable WA states are constant ^ lines. § = 95.5 and § = 660.7 correspond to the middle of the WA states of the Standard curve at 
10^ K and 10*' K respectively and g = 3630.8 corresponds to the middle of the 10^ K phase for the Z/Zq = 3 stability curve. 



following subsections we have studied the effect of different chem- 
ical compositions of the absorber on the stability curve. The illu- 
minating continuum is assumed to be the same as that in the Stan- 
dard set values (Table 1). We have generated the stability curves 
for (i) super-Solar abundances, (ii) sub-Solar abundances and (iii) 
other abundances where specific elements or group of elements 
have been removed from the solar metallicity. 

7.1 Super-Solar Abundance 

The left panel of Figure 7 shows stability curves for gas having 
super-Solar abundances as defined in Hamann & Ferland (1993). 
In addition to the 10'' K and the 10* K phases, all the curves also 
exhibit WA states. However the temperatures of the WA phases and 
their distribution in the phase space change with metallicity. The 
nature of stability exhibited can be separated into three different 
groups showing distinctly different trends. 

All the curves in the group with solar to three times solar ex- 
hibit stable gas at T ~ 10^ K, with the extent of the stable region 
decreasing as the metallicity is increased. ^5 = 95.5 for the Stati- 
dard curve, being the ionization parameter corresponding to the 
middle of the 10'' K WA phase. The range of ionization parame- 
ters for the lO'^ K gas is 50 < 5 < 182 for the Standard curve. 
The ^ = 95.5 drawn in the left panel of Figure 7 indicates that 
the 10^ K WA is pushed towards higher values of ^ as the metal- 



licity increases. Another stable states appears in all the curves at 
T ~ 10^ K, but this time with the increase in metallicity, the ex- 
tent of the stable region is increased, = 660.7 for the Standard 
curve, where is defined to be the ionization parameter corre- 
sponding to the middle of the 10^ K stable WA phase which shows 
the range 380 < £, < 1202 for the Standard curve. The solid black 
line across the 10® K phases of the different stability curves is the 
^ = 660.7 line. Thus, we see that with the increase of metallicity, 
^6 is decreased slightly and the 10® K stable phase moves to lower 
^/T and higher T. A very narrow (in the range of log (^/T)) phase 
is seen to appear at T ~ 10^ K for metallicity of 3 x solar. The 
detailed quantitative comparisons are made in Section 8. 

The curves in the next group with 4 < Z/Zq < 7, are very 
similar to each other, with a very narrow stable state at 10^ K, 
which is not in pressure equilibrium with the 10® K phase and the 
10^ K state becomes more prominent. As abundances are increased 
further (Z/Zq > 7) the stable state at 10^ K disappears, while the 
10® K and 10*^ K states remain quite wide and in pressure equi- 
librium with one another. The stability curves for the super-Solar 
abundances thus suggest that the possibility for the warm absorb- 
ing medium to exist in discrete phases in pressure equilibrium with 
each other increases with increasing metallicity. We return to these 
issues in Section 8 with quantitative details. 

The stability curves for the super-Solar metallicities with 
Z/Zq > 3 are seen to have gas at T < 6 x 10* K to be in pres- 



© 0000 RAS, MNRAS 000, 000-000 



12 Chakravorty et al. 




Figure 8. Stability curves for different chemical compositions with certain individual elements or certain groups of elements, as labeled, set to zero. "X-ray 
Group Off signifies the absence of the group of metals C,0, Ne, and Fe in the gas abundance, whereas "o-elements Only" represents an absorbing gas having 
only o-elements (Ne, Mg, Si, S, Ai', Ca, Ti) as heavy elements. See text for further details. The curve labeled 'Standard' in the two panels uses the Standard 
set of values (Table ??). We have also plotted the zero metals curve in both the panels for comparison 



sure equilibrium with the WA states at higher temperatures. Such 
stable states at ~ 10* K are also highhghted in the corresponding 
stabihty curves in the left panel of Figure 7. Gas at ~ 10* K might 
be representing the 'BELR clouds'. Hence this phenomenon might 
indicate the interesting possibility of some connection between the 
X-ray and UV absorber, as discussed by Turner etal. (1995) and 
Elvis (2000). However, these results have to be checked out further 
with more realistic AGN spectra which is beyond the scope of this 
paper and will be addressed in our next publication (Chakravorty et 
al. , in preparation). 

We have also investigated the possibility of high metallicity 
gas making it possible for WA to exist even when the X-ray ionizing 
continuum is very flat. The stability curves are shown in the small 
inset window of the left panel of Figure 7. Even for solar abun- 
dance, the a — 0.2 curves did not show any significant 10^ K WA 
and higher metallicity diminishes the chances further. At higher 
temperatures of ~ 10*' K, very narrow regions of stability appear 
for Z/Zq > 5. However, even in these cases there is no multi- 
phase exhibited by the stability curves. Thus even if the gas has 
high metallicity, WAs are not likely to be present in the intervening 
medium, if the ionizing continuum is very flat. 



7.2 Sub-solar abundance 

Gas in AGN has super-Solar metallicity (Hamann & Ferland, 
1999). However, AGN illuminate low density gas on large scales 
which has low abundances, especially at high redshift. Sub-Solar 
chemical composition are relevant for such line of sight absorbing 
systems including Damped Lyman-a systems in the inter galac- 



tic medium and smaller galaxies like the Magellanic clouds. The 
inter-cluster medium, for example, has ~ 0.3^© abundance and in 
certain systems there are even stars with an abundance of O.OIZ© 
(Pettini, 2006 and references therein). It is, therefore, interesting to 
examine the effects of sub-Solar abundance on the stability curve. 
We show the low metallicity stability curves in the right panel of 
Figure 7. As before, the solid black lines across the stable WA 
phases are the constant ^ lines. 

The ionization parameter ^5 corresponding to the middle of 
the 10^ K stable phase increases from its standard value of ^5 = 
95.5 with the decrease in metallicity and there is gradual enhance- 
ment of the thermal stability of gas. On the other hand, with almost 
constant ~ 660.7, the stable phase at log T ~ 6 moves towards 
lower temperature and higher pressure (^/T) and its extent in the 
log(^/r) shrinks until such a phase disappears for O.OZq. Detailed 
quantitative comparisons will be presented in Section 8. 

Absorbing gas with zero metals might be relevant for super 
massive black holes (SMBHs) or their seeds at redshift z > 6, 
which are formed by direct collapse and where no star formation 
precedes the radiating SMBH (Volonteri & Rees, 2005; Begelman, 
Volonteri & Rees, 2006). We have shown the zero metals stability 
curve in the right panel of Figure 7. There is a single phase WA at 
10® K with significantly high ^5 — 437 as compared to the Stan- 
dard ^5 = 95. 



7.3 Influence of individual elements and groups of elements 

Figure 6 has shown how the interaction between gas and radiation 
is sensitive to the ionization potentials of various ions of elements 
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like C, O, Fe and Ne indicating that thie WA states are likely to 
be affected not only by the overall abundance of the absorber, but 
also by specific elements and groups of elements. Figure 8 shows 
the effect on the stability curves if certain elements or groups of 
elements are not present in the absorbing gas. In both the panels of 
the figure, we have also drawn the stability curve for the Standard 
set of parameters and zero metals for comparison. 

The WA consists of absorption edges and lines of different 
ionization stages of oxygen, carbon, neon and iron. We define the 
X-ray group constituting of only these four elements. The left frame 
of Figure 8 shows the results of simulating various stability curves 
corresponding to the removal of different individual elements of the 
group from a solar composition gas. Setting carbon abundance to 
zero does affect the stability curve at logT < 4.5 and the devia- 
tions imply that Carbon acts as a cooling agent at those tempera- 
tures. However at higher temperatures relevant for WA, Carbon has 
no significant effect. The removal of oxygen and neon leads to de- 
viations in the stability curves which indicate that for T < 10^ K 
oxygen acts as a significant cooling agent with some contributions 
from neon, whereas for T > 10^ K both act as heating agents. 
Removal of iron has the maximum effect as it renders the curve 
relatively featureless in the range 4.5 < log T < 7. Such an ef- 
fect is expected because at higher energies iron is the most abun- 
dant absorbing element and plays a major role to stabilise the gas 
by acting as a significant heating agent. This result has interesting 
cosmological consequences. Type la supernovae created iron when 
the Universe was about 1 Gyr old (Hamann & Ferland, 1999 and 
the references therein). Hence there would be an interesting transi- 
tion in the nature of the WA at this age; WAs devoid of iron from 
before this era cannot show a significant 10® K phase. Future mis- 
sions like Constellation X might be able to detect such a change. 
Like the removal of iron, the removal of the X-ray group (O, C, Ne, 
Fe) as a whole zlso flattens out the stability curve, but to a larger ex- 
tent as shown in the right panel of Figure 8. If the absorbing system 
is devoid of the X-ray group, the stable states at 10^ K are pushed 
towards higher values of ^ with ^5 = 158 as compared to ^5 = 95 
for the Standard curve and there are no higher temperature stable 
states thus nullifying the possibility of multiphase WA. 

In addition to the elements of the X-ray group, it is instructive 
to study the effect of the Q-elements (Ne, Mg, Si, S, Ar, Ca, Ti) 
produced together with oxygen in type II supernovae which pre- 
cede the type la supernovae where iron and associated elements 
are formed. Thus the relative abundance of these elements gives a 
handle on the stellar evolution history and the age of the galaxy. 
Ne is the only a-element which has been studied as a part of the 
X-ray group (O, C, Ne, Fe) as well, and its effect on the stabil- 
ity curve has already been discussed above. Among the remain- 
ing a-elements, magnesium, sulphur, argon, calcium and titanium 
have negligible effects. Only silicon is a significant cooling agent at 
log T ~ 5, so that removal of silicon shows an almost vertical jump 
from 10® K to 10® K in the corresponding stability curve in the left 
panel of Figure 8. As a group however, a-elements suggest inter- 
esting consequences. According to the evolution history of metals 
in the Universe, the a-elements are formed first, followed by the 
other metals in subsequent stages. Hence we investigated the na- 
ture of an absorbing gas which has only the a-elements as metals 
and the corresponding stability curve is shown in the right panel of 
Figure 8 labeled as 'a-elements Only'. For log T > 5.2 the curve 
is same as the "X-ray group removed" curve and at lower tempera- 
tures, log T ~ 4.5, it traces the 0.0^© curve. These results suggest 
that a certain amount of evolution of metals is necessary for any 
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Figure 9. A summary plot depicting the stability curves for a range of chem- 
ical compositions. The curves marked 'Fe' and 'X-ray Group Off' corre- 
spond to absorbers where iron and the X-ray group elements respectively 
are absent. An absorber with a composition where only the o-elements are 
present as the metals, is represented by the stability curve labeled with 'a- 
elements On'. There is a gradual flattening of the stability curve as one 
moves from 5Zq to O.OZq. 



absorbing systems to show WA signatures and might set an upper 
limit to the redshifts at which warm absorbers can form. 



7.4 Summary of abundance effects on stability curves 

In Figure 9 we have collected the curves showing the major ef- 
fects due to varying the chemical composition of the absorber. As 
the metallicity of the absorber is decreased from 5Zq to Zq and 
to QAZq, the stability curve flattens gradually. The possibility of 
finding multiphase WA decreases with decrease in metallicity. This 
result is in conjunction with those found by Komossa and Mathur 
(2001). It also decreases if the absorbing medium is deficient in 
some elements. If a WA lacks iron, but has all the other elements 
with solar abundance, it then has nearly the same stability proper- 
ties as an absorber with half solar abundance. The X-ray group of 
elements are the most influential in making multiple phases pos- 
sible. Absorbing gas at early epochs of star formation when only 
a-elements formed show the same thermal properties as absorbers 
devoid of the X-ray group of elements, thus suggesting that the 
properties of WA are significantly dependent on the evolution his- 
tory of metals in the Universe. Figure 9 also shows the stability 
curve for O.OZq; such a primordial gas of hydrogen and helium 
has a single phase of absorbing gas at 10® K for a significantly 
high ionization parameter No stable state at 10® K and hence no 
multiphase solutions are possible if the metallicity of the absorber 
is < 0.2Zq. Metals thus seem to be necessary for absorbing gas to 
exist at 10® K < T < 10^ K. 

Hess et al. (1997) have carried out an extensive study trying 
to account for the unstable thermal equilibrium zones in tempera- 
ture ranges relevant for AGN and X-ray binaries paying attention to 
even individual ions of various metals responsible for causing the 
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instabilities. They have predicted the ionization structure and ther- 
mal states of absorbers with a wide variety of chemical composi- 
tions. The results posed in this paper especially in Section 7.3 agree 
qualitatively with their results but as expected, they vary in quanti- 
tative details. For example, like them, we also find oxygen and iron 
to be playing crucial roles in the ionization and thermal balance 
equations. However, unlike Hess et al. in this paper we find that 
oxygen and iron influence distinctly different temperature zones of 
the stability curve. As mentioned in their paper, such minute details 
would crucially depend upon the accuracy of the atomic physics in- 
corporated which is definitely addressed much better in the current 
photoionization codes than in 1997. 

The determination of the particular ions responsible for the 
heating or cooling of the gas at WA temperatures, is very sensitive 
to the exact shape of the ionizing continuum between 50 to 500 
eV. Such investigations and quantitative comparisons are beyond 
the scope of this paper and will be addressed in our subsequent 
publications (Chakravorty et al. , in preparation). 



8 MULTIPHASE ANALYSIS 

The equilibrium curves we have derived in the previous sections 
often contain segments that allow phases at different temperature 
to occur at similar pressure ^/T. This has clear physical interest as 
distinct phases in pressure equilibrium have been derived in a num- 
ber of cases of WAs (Krongold etal. , 2003, 2007, Netzer etal. , 
2003), while in others it has been claimed that there is a continu- 
ous range of ionization parameter (Ogle etal. , 2004, Steenbrugge 
etal. , 2005). In this section we quantify the extent of these mul- 
tiphase regions. However, we would like to note that for carrying 
out such an exercise, the atomic database needs to be accurate. As 
mentioned a number of times in the previous sections, care has been 
taken to upgrade the tmderlying atomic physics in cloudy with the 
most recent updates. 

In Table ?? we present detailed quantitative estimates of the 
WA properties for the various parameter sets used in this paper. 
For convenience of representation we introduce a few definitions 
for use in Table ??. The first four columns describe the parameter 
set used to derive the stability curve and the WA properties. ^„ in 
columns 5, 9 and 13 and r„ in columns 6, 10 and 14 respectively 
are the ionization parameter and the temperature corresponding to 
the middle of the 10" K stable WA phase, where T„ is expressed 
in units of 10" K. A„[log(^/T)] in columns 7, 11 and 15 give 
the range of log(^/r) over which stable states of WA are exhib- 
ited for a temperature of 10" K. For WA to exist in discrete multi- 
ple phases, different temperature stable states require to exhibit the 
same values of^/T. The larger the range of these common values, 
higher we presume to be the probability of a multiphase WA. Such 
intersection in the range of ^/T is represented by Ase [log(^/T)] 
for 10^ K and 10® K phases and by A67[log(C/r)] for the states 
at 10^ K and 10® K and their values are noted in columns 8 and 12 
respectively. Np in the second last column is the number of stable 
states in the stability curve in the temperature range 5 < log T < 7 
and the number of states in pressure equilibrium with each other is 
given by Nmp in the last column. We would note that the values 
of the various quantities recorded in Table ?? are useful for un- 
derstanding trends in the variation of the physical quantities as the 
parameter sets are changed. However, for more robust estimations 
of these nvunbers and specific predictions of the WA properties, we 
require more detailed modeling of the ionizing continuimi and the 



gas abundances. Such investigations will be described in a subse- 
quent publication (Chakravorty et al. in preparation). 

In Figure 10, we plot A„[logT] against A„[log(f/r)] for 
a range of continuum shapes (left) and abundances (light), where 
A„[logT] and A„[log(C/T)] are the range of T and ^/T for the 
stable states at 10" K. Presumably, the higher the values of these 
two quantities, the better are the chances of existence of the state 
at 10" K, because any perturbed gas is more likely to find a larger 
stable region. The various parameter sets used to describe WA con- 
ditions in this paper have been represented by different symbols 
which have been explained in the bottom panel. The gray square 
represents the Standard set of parameters and the rest of the sym- 
bols indicate the changes from it. 

From the systematic analysis we can see the following trends: 

• The variation of soft X-ray slope a of the ionizing continuum 
shows that WAs exist only if a > 0.2 and will have higher temper- 
ature gas if a > 0.5. 

• Increases in the soft X-ray slope a decrease the ionization pa- 
rameter of the 10'' K WA while decreasing its temperature, but in- 
creases both 5 and T for the 10® K phase. Both A5[log(^/T)] and 
A6[log(^/r)] increases with the value of a, which is also demon- 
strated by the square data points in the top and the middle left pan- 
els of Figure 10 accompanied with increase in A„ [log T)], indicat- 
ing that WA phases stabilise with steeper soft X-ray continua. 

• The Multiphase scenario, however, is not faciUtated sig- 
nificantly with the variation in a. Multiphase solutions at 
10^ and 10® K are seen only for the standard value of a = 0.8. 
Although the a — 1.1 stability curve also shows 2 phases, the 
value of A56[log(^/T)] is very low, making the possibility of mul- 
tiple phases in equilibrium seem unlikely. 

• We have found very interesting difference between our results 
and those obtained by Krolik & Kriss (2001). They have drawn 
stability curves for ionizing continuum having properties similar to 
our Standard set SED having Oox ~ 1.2 and a = 0.8. The stability 
curves shown in KKOl have a near vertical jump in temperature 
from ~ 3 X lO"* to ~ 10® K as ^ runs from ~ 40 to ~ 1500 thus 
exhibiting a continuous distribution of stable states in ^/T — T 
space. For the same set of parameters (ie, the Standard set) we not 
only get two distinct WA phases in pressure equilibrium with each 
other, but also do not see any drastic change in temperature over 
a narrow range of log(^/T), although our ^ values also run from 
~ 50 to ~ 1200. 

• As the value of the UV-X-ray slope aox is increased, the ion- 
ization parameter for the 10^ K WA increases steadily accompa- 
nied by a small decrease in temperature, also increases but with 
a small increase in temperature. 

• The range A5[log(^/r)] remains almost constant for 1 < 
ctox < 1.6, but abruptly increases for the value of 1.8. 
A6[log(^/T)] grows abruptly for Oox ~ 1-6 before which it in- 
creases steadily but slowly with the increase in aox- There is no 
significant WA state at 10® K for aox = 1-8, Comptonization be- 
comes the dominant heating and cooUng processes at these temper- 
atures for such an SED. The triangles in the left panels of Figure 10 
show these trends. 

• Multiphase solutions at 10^ and 10® K are seen in the range 
1 < Oox < 1.6. For aox = 1-4, there is an additional stable state 
atlO'^ K,but A67[log(C/r)] is very low to claim any possibiUty of 
this state being in pressure equilibrium with the lower temperature 
phases. 

• Through the increase of metallicity of the absorbing gas from 
0.2^0 to 9Zq, the 10* K WA maintains almost constant tempera- 
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Table 2 The physical and thermal stability properties of WAs corresponding to various parameter sets considered in this paper. The various physical physical properties for the WA are given 
in coluirms 1 - 4. Colunms 5-7 note the presence of the different components of the WAs at T ~ 10^, 10®, and lO'^. The number of WA phases is given by colunm 8 while a qualitative 
discussion of the WA is given in column 9. 
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Figure 10. For the stable states of the WAs we have plotted A[logT] against A[log(^/T)] and have considered the parameter sets where we have varied 
the SED of the ionizing continuum and the metallicity of the absorbing gas. See text for details. Various symbols used to identify the different parameter sets 
used in this paper, are explained in the bottom panel. The gray square represents the Standard parameter set and for the others the deviations from that set are 
mentioned. 



ture. For sub-Solar metallicity, ^5 decreases steadily through the 
increase of abundance, but this trend is inverted for super-Solar 
metallicity. The 10*' K WA, on the other hand, shows a rise in tem- 
perature Te throughout the metallicity variation, while remains 
constant within a factor of 1.4. 

• There is an anti-correlation between A5[log(^/r)] and 
A6[log(^/r)]; the 10^ K WA shrinks in log(^/r) until it van- 
ishes for QZq, whereas the phase at lO*" K progressively grows in 
log(^/r) from OAZq to 9Zq. This trend is also depicted in the 
top and the middle right panels of Figure 10 by the open stars and 
the closed circles. A lO'^ K phase starts appearing from 3Zq , and 
A7[log(^/r)] increases with abundance. 

• Interesting multiphase properties are exhibited by WA with 
varying abundance. A56[log(5/r)] is significant in the range 
0.4 < Z/Zq < 3 for which lO'^ and lO" K gas are in pressure 
equilibrium with each other. Beyond this range for Z/Zq > 5, 



the stability curves show the interesting feature of significant mul- 
tiphase solution between the 10® and 10^ K phases; a property de- 
picted by only these three stability curves among all the parameter 
sets considered in this paper. 



• In the case of the parameter sets having super-Solar metal- 
licity, we find an interesting phenomenon depicted by the stability 
curves. For Z/Zq > 3, the WA stable phases (5 < logT < 7) 
are seen to be in pressure equilibrium with gas at ~ 10* K, which 
might be representing the 'broad emission line region (BELR) 
clouds' . Such a feature in the stability curve might indicate the con- 
nection between the X-ray and UV absorber. 
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9 SUMMARY 

Absorption edges due to various partially ionized elements are seen 
in the soft X-ray spectra of about half of the Seyfertl galaxies and 
observed quasars. These absorption features are due to the pres- 
ence of a warm (~ 10^ — lO'''^ K) gas along our line of sight to 
the AGN. We have shown that the existence and nature of this ab- 
sorption component crucially depends on the shape of the ionizing 
continuum from the central source which illvuninates the gas, and 
also on the chemical composition of the absorbing medium. We 
have used the stability curve as a tool to study the stabiUty of the 
WA as a function of its various physical properties. 

We have considered an ionizing continuimi with two power- 
law components which together reproduce the observed X-ray 
spectral index as well as the observed values of aox- Systematic in- 
vestigation shows that if the X-ray spectral index, a, of the ionizing 
continuum intervening gas of solar metallicity is flat with a ~ 0.2, 
then the AGN cannot have a WA. For moderately steep X-ray spec- 
tra withO. 5 < a < 1.1, it is possible to sustain a WA. A multiphase 
nature of the WA is, however, seen only for a ~ 0.8. It is an in- 
teresting coincidence that most of the observed quasars also have 
soft X-ray slopes similar to a ~ 0.8 and in future studies we shall 
investigate whether such a coincidence leads to deeper physical im- 
plications. If the ionization continuum becomes even steeper then, 
instead of discrete phases, the absorbing gas exhibits a continuous 
distribution of ^/T and T. We have also verified that the nature of 
the WA is not affected by the higher energy cut-off of the ionizing 
continuum. 

Temperature estimates for warm absorbing gas were generally 
accepted to be independent of the density of the gas, because the 
ionization and thermal equilibrium is believed to be achieved by 
the balance between photoionization and recombination. We show 
in this paper, in agreement with Rozanska et al. (2008), that the 
stability properties of highly dense WA might show density depen- 
dence influenced by sufficiently soft UV dominated ionizing con- 
tinuvun. This is a very interesting result because it implies that for 
such objects, we would not require independent estimates of den- 
sity from observations; modeling the UV spectra itself can provide 
constraints on gas density. 

The chemical composition of the absorbing medium plays a 
critical role in determining its multiphase nature. The higher the 
metallicity of the medium, the greater is the probability of a mul- 
tiphase parameter space. The role of individual elements or cer- 
tain important groups of elements were also studied. Amongst the 
elements, iron, which is formed from type la supernovae no ear- 
lier than 1 Gyr, plays the most important role as a heating agent 
at higher temperatures (> 10^ K) while oxygen is a significant 
cooling agent for T < 10^ K. The X-ray group (C, Ne, O, Fe), 
tmderstandably, has significant influences, because the constituent 
elements are the ones which have important atomic transitions in 
the energy range relevant for warm absorbers (0.3 — 1.5 keV). If 
the absorber is over abundant in iron, oxygen or in the X-ray group 
of elements then the absorbing medium is likely to be multiphase. 
A very interesting result seems to suggest itself from the study of 
influence of a-elements (Ne, Mg, Si, S, Ar, Ca, Ti), which are the 
first metals formed in the Universe along with oxygen through the 
supernovae of type II. The thermal properties of the warm absorber 
change drastically from the time when the absorbing gas is essen- 
tially constituted of only a-elements as metals to the time when 
the absorber is enriched enough to have the whole X-ray group 
included in its chemical composition. We have also observed that 
zero metal abundance gas does not produce WAs in AGN. Our anal- 



ysis on the chemical composition of the absorbing gas shows that 
the warm absorber temperature range is very sensitive to the exact 
abundance of the absorber and there is an indication that the warm 
absorber might have super-Solar metallicity. While fitting soft X- 
ray data and determining the physical properties of the warm ab- 
sorber in individual objects care should be taken to include the ef- 
fect of super-Solar metallicities. 

When we consider super-Solar metallicity, the stability curve 
shows a 10^ K stable state is in pressure equilibrium with the WA 
stable phases. This might imply some connection between the WAs 
and the BELR. We plan to take up this issue in future work. 

We have used an upgraded version of cloudy (C07.02) where 
the most recent developments in atomic physics have been incor- 
porated. This gives us the confidence to carry out a quantitative 
analysis on the multiphase nature of the WA presented in rigorous 
details in Section 8. We have been able to provide specific trends 
of the behaviour of the WA as a function of the prevalent physical 
conditions like the shape of the ionizing continuvun and the chem- 
ical composition of the gas. Comparison between observations of 
WA properties and simulations can be made more concrete by pro- 
viding such quantitative measures. We will extend this exercise to 
greater details in our subsequent pubUcations. 

In this paper we have used simple broken power-laws to define 
the ionizing continuum for the absorbing gas. However, the soft 
excess in Ultraviolet is often modeled using blackbody components 
which peak at T ~ 150 eV. Moreover the continuum spectra of 
AGN also has the disk blackbody component at 20 eV < T < 
30 eV. We intend to investigate the effects of such continua on the 
WA properties as the next step in this systematic analysis. 
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